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Abstract 

We present the exact solution of two-body motion in (1+1) dimensional dilaton 
gravity by solving the constraint equations in the canonical formalism. The determining 
equation of the Hamiltonian is derived in a transcendental form and the Hamiltonian is 
expressed for the system of two identical particles in terms of the Lambert W function. 
The W function has two real branches which join smoothly onto each other and the 
Hamiltonian on the principal branch reduces to the Newtonian limit for small coupling 
constant. On the other branch the Hamiltonian yields a new set of motions which can 
not be understood as relativistically correcting the Newtonian motion. The explicit 
trajectory in the phase space (r,p) is illustrated for various values of the energy. The 
analysis is extended to the case of unequal masses. The full expression of metric tensor 
is given and the consistency between the solution of the metric and the equations of 
motion is rigorously proved. 
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1 Introduction 



One of the oldest and most notoriously vexing problems in gravitational theory is that of de- 
termining the (self-consistent) motion of N bodies and the resultant metric they collectively 
produce under their mutual gravitational influence [[TJ. In general there is no exact solution 
to this problem (although approximation techniques exist 0) except in the case N = 2 for 
Newtonian gravity, since energy dissipation in the form of gravitational radiation obstructs 
attempts to obtain an exact N = 2 solution. Only the static sector of the Hamiltonian has 
thus far been determined exactly ||. 

Lower dimensional theories of gravity do not contain gravitational radiation and so offer 
the possibility of making useful progress on this problem. For example in (2 + 1) dimensions, 
the absence of a static gravitational potential allows one to generalize the static 2-body 
metric to that of two bodies moving with any speed [|J. In (1 + 1) dimensions one must 
necessarily consider a dilatonic theory of gravity [f| since the Einstein tensor is topologically 
trivial in two dimensions. One such theory in this class has been of particular interest insofar 
as it has a consistent Newtonian limit |J (a problematic issue for a generic dilaton gravity 
theory [[?]]), allowing for the formulation of a general framework for deriving a Hamiltonian 
for a system of particles |§. In the slow motion, weak field limit this Hamiltonian coincides 
with that of Newtonian gravity in (1 + 1) dimensions. 

Motivated by the above, we consider here the problem of the relativistic motion of 2 
point masses under gravity in (1 + 1) dimensions. We work in the context of the dilatonic 
theory of gravity mentioned above 0. Both the classical and quantum properties of this 
theory (referred to as R = T theory) have been extensively investigated || ^ [II], |TTJ , and it 



contains the Jackiw-Teitelboim lineal gravity theory |12| as a special case ||. The specific 
form of the coupling of the dilaton field \& to gravity is chosen so that it decouples from the 
classical field equations in such a way as to ensure that the evolution of the gravitational field 
is determined only by the matter stress-energy (and reciprocally) || ||. It thereby captures 
the essence of classical general relativity (as opposed to classical scalar-tensor theories) in 
two spacetime dimensions, and has (1 + l)-dimensional analogs of many of its properties 
H [II]]. Furthermore, this theory can be understood as the D — > 2 limit of general relativity 
(as opposed to some particular solution(s)) |T3| . 

We consequently find that the motion of the bodies is governed entirely by their mutual 
gravitational influence, and that the spacetime metric is likewise fully determined by their 
stress-energy || [|. Unlike the (2 + 1) dimensional case, a Newtonian limit exists, and there 
is a static gravitational potential. The solution we obtain gives the exact Hamiltonian to 
infinite order in the gravitational coupling constant. Hence the full structure of the theory 
from the weak field to the strong field limits can be studied. While some of the phase-space 
trajectories we obtain can be viewed as relativistic extensions of Newtonian motion, we find 
that for sufficiently large values of the total energy a qualitatively new set of trajectories 
arises that cannot be viewed in this way. 

The outline of our paper is as follows. In section 2 we recapitulate the derivation of the 
canonically reduced iV-body Hamiltonian in (1 + 1) dimensions. In sections 3 and 4 we solve 
the constraint equation and then derive an expression for the exact Hamiltonian in the 2 
body case. In section 5 we analyze the motion in the case of equal masses, and in section 
6 we consider the unequal mass case. In section 7 we solve for the spacetime metric and in 
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section 8 we investigate the test-particle limit of our solution. In section 9 we consider the 
the dependence of the Newtonian limit on dimensionality. We close our manuscript with 
some concluding remarks and directions for further work. 



2 The canonically reduced TV-body Hamiltonian 

First we shall review the outline of the canonical reduction of (1 + 1) dimensional dilaton 
gravity ||. The action integral for the gravitational field coupled with N point particles is 



a=l J I 



dzt dz h „ 



1/2 



dr n . dr n 



[x 



Za{r a )) 



where \I/ is the dilaton field, g^ is the metric (with determinant g), R is the Ricci scalar 
and r a is the proper time of a-th particle, with k = 8ttG/c 4 . The symbol V M denotes the 
covariant derivative associated with g^ v . 

The field equations derived from the action (JJ) are 



R - g^V^V^ = 



d j . .dz v a \ 1 



dr. 



dr„ 



dz^dz^ 
dr a dr a 



. 



where the stress-energy due to the point masses is 



N 1 

a=l J V 9 



(j ® vp dr n . dr, 



d7 a d?P 

'-" (l ' a 5 2 (x-z a (r a )) . 



Inserting the trace of @ into (§) yields 



R = kT» 



(2) 
(3) 
(4) 

(5) 
(6) 



Particle dynamics in R = T theory may therefore be described in terms of the equations 
and (H), which forms a closed system of equations for the gravity /matter system. 

At first sight it may seem that the dynamics is independent of the dilaton field, since 
both (H) and @ do not include Note, however, that all three components of the metric 
tensor cannot be determined from (^j), since it is only one equation. The two extra degrees of 
freedom are related to the choice of coordinates. If the coordinate conditions are chosen to be 
independent of equation (|5|) determines the metric tensor completely. However this need 
not be the case, and so more generally we need to know the dilaton field through which 
the metric tensor is (indirectly) determined - it is this field that guarantees conservation of 
the stress-energy tensor via (|3]). 
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In the canonical formalism the action (|T|) is written in the form 

1 = Jdx 2 !^p a z a 5(x - z a (x )) + 71-7 + + N R° + N^ 1 



(7) 



where 7 = gn, No = (— g 00 ) 1 / 2 ,Ni = gio, it and IT are conjugate momenta to 7 and $ 
respectively, and 



^y/l « Vv^/ a V 7 

(8) 

i? 1 = -TT--U^' + 27v' + y^5(x-z a (x )) . (9) 

7 7 a 7 

with the symbols ( ' ) and ( ' ) denoting <9 and di, respectively. 

The transformation from (|]) to (0) is carried out by using the decomposition of the scalar 
curvature in terms of the extrinsic curvature K via 

^g-R = -2d (^K) + 2d 1 [ v ^(N 1 K - ^d^)} 

where K = (2No'~f)~ l (2diNi — r y~ l Nid\ r ) — <9o7) and changing the particle Lagrangian into 
first order form. 

The action (^) leads to the system of field equations : 



7T 



+ n x + - + Y,-2^~ *»(* ))) + Kt^*' + Ni- = o 

{ T 7 a 7 J 2/c^7 7 7 

(10) 

V 

7 - iV (2«*/777r - 2«V7n) + ^i— - 2iV( = (11) 

7 

R° = (12) 
i? 1 = (13) 

n + ^(-I^n + — J-pW + -^—Nfi = (14) 

7 2k^7 k^7 

* + JVb(2« A /77r)-^i(-*') = (15) 
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4-iVo / + — = • (17) 

ry 1 < < 
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In the equations (fL6l) and (|T7|), all metric components (No, N%, 7) are evaluated at the point 

df(x) 



x = z a and 



df_ 

dz a 



dx 



This system of equations is equivalent to the set of equations (||), (|3|) and (|J). 

The action ([?[) also shows that iVo and Ni are Lagrange multipliers, and the equations 



(|i~2"|) and (|T^) are constraints. We may investigate the canonical structure of the theory via 
the generator which arises from the variation of the action at the boundaries : 



G = J dx ^2p a 3(x — z a)^ z a + ffSh — ^SH 



;is) 



where h = 1 + 7. This form was obtained by adding a total time derivative — do(U^>) to 
the original action and taking the constraints into account. Since the only linear terms 
in the constraints are fif'/y/j)' and it', we may solve for these quantities in terms of the 
dynamical and gauge (i.e. co-ordinate) degrees of freedom. Bearing this fact in mind, we 
transform the generator fll8|) to 



G = J dx l^2p a 6(x - z a )5z a 



1 



2tt' 



h-l 



K \yfh 



s 



K 

A 



h-l 



-IN' 



Vh 



1 



S {2A h '. 



(19) 



where 1/ A is the inverse of the operator A 
and we have discarded surface terms. 
Adopting the coordinate conditions 



d 2 /dx 2 with appropriate boundary condition, 



x 



4^' 

2A 



t 



1 



allows the generator ([19]) to be expressed in the canonical form 



G = J dx jx)Pa<H 



where 



T 



x - z a )5z a 



H = — 
2n' . 



(20) 

(21) 

(22) 
(23) 



In deriving this expression for T 0/1 , we made use of the differential forms of the coordinate 
conditions 

7 = 1 n = . (24) 
As discussed in M, consistency between the integral form (B^) of the coordinate conditions 



and the differential form (p4|) is assured when one retains the appropriate boundary condi- 
tions for the integral operator 1/ A and takes a limiting procedure by introducing a regulator. 
Thus, TC is the Hamiltonian density of the system and T 01 is the momentum density. 
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The action (^) is similarly tranformed and reduces to 



J dx 2 V£2p a z a 8(x - z a ) . 



(25) 



Thus the reduced Hamiltonian for the system of particles is 



H = I dxH = -- I dxA& (26) 



K 



where ^ is a function of z a and p a and is determined by solving the constraints which are 
under the coordinate conditions (p4|) 

1 



A$ - -(^') 2 + K 2 vr 2 + \jvl + ml5(x -z a )=0 (27) 

^ a 

2^ + = . (28) 

a 

The expression of the Hamiltonian ( pB"D is analogous to the reduced Hamiltonian in (3 + 1) 
dimensional general relativity. The proof of the consistency of this canonical reduction was 
given in || : namely the canonical equations of motion derived from the reduced Hamiltonian 
(HH) are identical with the equations and fll7j) . 



3 Matching conditions and the solution to the con- 
straint equations 

The standard approach for investigating the dynamics of a system of particles is to derive 
an explicit expression for the Hamiltonian, in which all information on the motion of the 
particles is included. In this section we solve the constraints (|27| ) and ( p8|) for the system of 
two particles and determine the Hamiltonian (|28|). 
Defining and x by 

tf = -41og|0| n = X ' (29) 



the constraints (P7|) and (|28|) for a two-particle system become 



A0-^-( X ') 2 = j + ™? - + VpI + ml cj>{z 2 )8{x - z 2 ) } (30) 



4 uw r 4 

Ax = -- - z x ) + p 2 S(x - z 2 )} . (31) 



The general solution to fl31~|) is 



X = —g {pi | x - Z\ | +p 2 | a; - z 2 |} - eXrc + eC x . (32) 

The factor e (e 2 = 1) has been introduced in the constants X and C x so that the T-inversion 
(time reversal) properties of x are explicitly manifest. By definition, e changes sign under 
time reversal and so, therefore, does x- 
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Consider first the case z 2 < Z\, for which we may divide spacetime into three regions: 

Z\ < x (+) region 

z 2 < x < Z\ (0) region 
x < z 2 (-) region 



In each region \' is constant: 



X 



-eX-\{ Pl +p 2 ) 
-eX + \( Pl -p 2 ) 
-eX + \{ Pl +p 2 ) 



(+) region 
(0) region 
(-) region 



(33) 



General solutions to the homogeneous equation A0 — (k 2 /4)(x') 2 <p = in each region are 

f </> + (x) = A + e^ x+ i {pi+P2) } x + B + e-^ x+ ^ 1+p ^ x 
(f> {x) = A e^{ x -^ Pl - p2) } x + B e-^{ x -i {pi ~ P2) } x • (34) 
6Jx) = A_e%{ x -^ {pi+P2) } x + B^e~^{ x --± {P1+P2) } x 



For these solutions to be the actual solutions to Eq.(|30D with delta function source terms, 
they must satisfy the following matching conditions at x = Z\, z 2 : 



(f)+(zi) = (f)o(zi) = <f>(zi) 
4>-{z 2 ) = (j) (z 2 ) = (j)(z 2 ) 

4>'+izi) - 0o(^i) = -\jv\ + rn\(j){zx) 

<P'q{z 2 ) - <p'_(z 2 ) = 

The conditions ( |35a|) and (|35c|) lead to 

gf { X+ % (Pl+P2)}zi^ + _|_ g-f {^+f (Pl+P2)}^l _g + 



'p 2 2 + m 2 2 <j){z 2 ) . 



(35a) 
(35b) 

(35c) 
(35d) 



(36) 



and 



,f {x+l(pi+ P2 )}z 1A 



f {^+f(Pl+P2)}2 



pf + + 2 {X - f ( Pl -P 2 )} ^ {X _, (P1 _ P2)U 
2{x + |( Pl +p 2 )} 

» + m\ - 2 {X - - p 2 )} r ., {x _ 4(pi _ p2)}zi 
2{x + l( Pl +p 2 



Bn 



(37) 



and then 
A, 

B 4 



p\ + mf + AX + ep 2 _«£ 



4X + e( Pl + p 2 ) 
4X + e( Pl +p 2 ) 



A + 



Pi+ml + epi 



AX + e( Pl +p 2 ) 



e - K (x+f P2 > 15o 



21 



An 



+ m\ - 4X - ep 2 « 
4X + e{ Pl +p 2 ) 



P1Z1 B 



(38a) 
(38b) 
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Similarly the conditions ( [35b|) and ( |35d|) lead to 



and 



e f { x ~i (pi+P2)}z2a__ _|_ e -f {x-f (pi+ P2 )}^ j g_ 
= et{ x -f^-P 2 )}^Ao + e-t{ x -f( pi -P 2 )} Z2 fi (39) 

g + mj - 2 {X - - {x _ |(pi _ p2)}z 
2{X-l( Pl +p 2 )} 



+ m 2 2 + 2{x - j(pi-p 2 )\ Kiy e, ,\ 

L ill lllip-f {*-f (Pl-P2)}*2 D 

2{X-f( Pl +p 2 )} 



(40) 



and then 



A_ = - V" 2 + ^ - 4X + £P1 e*~Aa - ^ + T 2 + ^ e-<*-^h Bo (41 
AX-e( Pl +p 2 ) AX-e( Pl +p 2 ) 



a) 

'jh+Pz) 

B- = v^ + ^-^ eK x- fP 0^ o+ V^ + ^ + 4X-e Pl e „ fp2Z2 
4X-e(pi+p 2 ) 4X-e( Pl +p 2 ) 

Since the magnitudes of both and x increase with increasing \x\, it is necessary to impose 
a boundary condition which guarantees that the surface terms which arise in transforming 
the action vanish and simultaneously preserves the fmiteness of the Hamiltonian. 

In the iterative analysis in || this condition has been shown to be 

fy 2 — 4k 2 % 2 = in the region Z\ < x and x < z 2 . (42) 

Since 

x _ I ~\ eX + l(Pi+P2)} x + eC x + \(Pi z i+P2 z 2) (+)region 
V 1 ~ \ eX ~ + V2 n x + eC * ~ K^i^i + ^2) (-)region 
the boundary condition implies 



A_ = B + = (44) 

KE , _ KE , K 

-( 

2 



-logv4+ - — (ptzi + p 2 z 2 ) = logB- + —{pizx +p 2 z 2 ) = -C x . (45) 



The condition flHT) leads to 



A = V^ + mj + epa r _4 x _. (pi _ pa) > <a (4g) 



7 



and 



/pi + mf — epi 

From (E6|) and (]47f) we have 



= Vg! + mj-4X- ep 2 c _ K{x _. (pi _ p2)h 

B . L2 _l ~2 



lp{ + m\- ep 2 - 4X)(\/pl + m\ + ep x - AX) 

= (\/pi + ml- t Pl )^p 2 2 + ml + ep 2 ) e 4*-ito-i*)}(* 

On the other hand the condition (HoT) leads to 



rP2 21 



4X + e(pi+p 2 ) 4X + e(pi+p 2 ) 



'4X-e( Pl +p 2 ) 4X-e( Pl +p 2 ) 



Using the notation 



we obtain 



A 
B 



Mi 
M 2 

Li 
L 2 
L+ 
U 



Pi 



mf — epx 



pl + ml + ep 2 



AX 



m - \lvl + m l 



AX + ep 2 - 
AX + e( Pl +p 2 ) 
AX - epi + ep 2 
AX-e( Pl +p 2 ) 



Pi + ml 



(M 2 L 2 ) 1 / 2 



L 

Mo 



r \ 1/2 r 
-^0 



Mi 



from 



(0) and (H). Substituting (0) and fl52|) into fl38a|) and (111 

A + = 



we 



Ll^ 


1/2 


) e-t L +^ 


Mi, 




L 2 ^ 


1/2 K 

e s L - 22 


M 2/ 
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and then 




In the case of z\ < z 2 we have to interchange the suffices 1 and 2. The equation fl48|) 
which determines X then generalizes to 

(yjpl + mj- ep 2 - 4JC)(yJp% + ~m\ + ep x - AX) 

= (\jpl + ml- m){4v\ + m 2 2 + ep^e^'^ 1 ^}^ (58) 
where r = z\ — z 2 and p a = p a sgn(zi — z 2 ). 



4 Determining equation for the Hamiltonian and the 
canonical equations of motion 

Since the solutions of give 

m-*{x + fr + *>} and f£ "-{X-L^} <«, 



the Hamiltonian H is 



H = -- I (LrAty 

K 



k [ J -°° k l->oo\0 + (L) <f>-(-L)j 
AX (60) 



and so (153) becomes 



p\ + m\ - ep 2 - H){\Jpl + ml + epi - H) 

= (\fpJ+^ 1 -ep 1 )(^pJ+^ 2 + ep 2 )e^ {H - e{pl - p2)} ^ . (61) 

Eq. (|6T| ) is the determining equation for the Hamiltonian, whose solution yields H as a 
function of (pi,p 2 ,r). 

Expanding ( pi] ) in powers of k yields the perturbative solution 



H = \JPi + ml + \l p 2 + ml + -(^ p\ + m\^J p 2 2 + m 2 2 - p^) \r\ 



+ -£■( VP? + m i Vi - Viyvl + m 2 2 )(zi - z 2 ) 



{(\fPl + m l\/P2 + m 2 ~ 2piP2)(\/p? + m l + ^Pl+mf) 



p{ + m\pl + p\Jpl + m z 2 \ 



n 2 e 



^ |-2(y / pf + m\\jp\ + m 2 2 - Pxp 2 ){pi - P2) + m\p 2 - Pimj^ \r\(z 1 -z 2 ) 

(62) 

up to 0(k 2 ). This is identical with the Hamiltonian derived in the iterative method in ||. 
For the case of z 2 < z 1 Eg. (|6l"D is 



L X L 2 = M x M 2 ei UAzi - Z2) . (63) 
Differentiating Eq. (|63|) with respect to z\ leads to 

^ + L 2 ) = - Ll L 2 [r- + L^ . 

Then we have the canonical equation 

dH 



Pi 



dz\ 
~^L LiL 2 

L\ + L 2 — ^-L\L 2 



(64) 



and similarly 

P2 = i L oL 1 L 2 

Li + L2 — -J-L1L2 

Differentiating Eq.(|63|) with respect to p\ leads to 

'dH Pl \ r (dH \ T f e «r /<9# N 



We have also the canonical equation 

dH 
opi 



e ^ L o L i 1 ^ 

L1+L2- fLxL2 J p j _|_ m 2 



and similarly 



Z2 = _e + "T~rT — ^tj — / 2 2 • ( 67 ) 



It is evident that the Hamiltonian and the total momentum P = p\ + p 2 are constants of 
motion: 

H = P=p 1 +p 2 = . (68) 
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5 Hamiltonian of two identical particles 

In this section we shall try to solve Eq.flBT]) for a system of two identical particles. We may 
choose the center of inertia frame with p\ = —p 2 = p. Then Eq. (|6lD becomes 



(H-^/p 2 + m 2 -ep) 2 = ^p 2 + m 2 -ep) 2 e^ H - 2eii)lrl . (69) 

After setting 



eq. ( |69"D becomes 



H = dp 2 + m 2 + ep ri Y (P, r ) ( 70 ) 

v K\r\ 

Y 2 e 2Y = Z 2 e 2Z (71) 

where Z = ^(a/p 2 + m ' 2 ~ e P)- 

Equation ( |7T|) has three solutions shown in Fig.l. The trivial solution, Y = Z, yields 
H = 2ep, which is unphysical because it has no interaction term. The second solution (curve 
A-B-O-C-D) is represented by 

Y = W(-Ze z ) Z <z = W-\e- 1 ) (72) 

and the third solution (curve E-F-G) is represented by 

Y = W(Ze z ) Z < (73) 

where W(x) is the Lambert VT-function defined via 

y . e y = x y = W{x) . (74) 



Fig.l 

Solutions to The points B and C represent the extremal Z and Y values of 

W~ 1 (l/e) = 0.278 on the principal branch. 

In general W is complex and multivalued. When x is real, the function has two real 
branches shown in Fig.2. The branch satisfying —1 < W(x) ( a solid line ) is denoted by 
Wq(x) and is referred to as the principal branch. The branch satisfying W(x) < — 1 ( a 
broken line ) is denoted by W-i(x), and is real-valued only for — 1/e < x < 0. Then for 
— 1/e < x < the function is double valued. The principal branch is analytic at x = and 
has a derivative singularity at x = —1/e beyond which W(x) becomes complex. The series 
expansion of the principal branch is given by 

W (x) = £ { -^—x n . (75) 



71=1 



Fig.2 
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The Lambert W function 



The correspondence between the curves in Fig.l and those in Fig. 2 is 

curve A-B-O-C-D curve P-Q-O-R-0 

curve E-F-G -<=>- curve O-Q-P 

Since the physical domain of Z is Z > 0, the only physical solution is ([72]), which yields the 
Hamiltonian 



W (— | (\r\ y/p 2 + m 2 - eprj 


exp 




\Jp 2 + m 2 — eprj 


]) 


h 


\r 





H = \ p 2 + m 2 + epsgn(r) — 

K\r\ 

(76) 

This Hamiltonian is exact to infinite order in the gravitational coupling constant. We can 
thus view the whole structure of the theory from the weak field to the strong field limits. 

The weak field expansion has already been given in the general case in flB2"p. The small 
p expansion is 



ran 


r 




-Sw 


K 


r 





ran 


r 


+ 2mK,\r\w + 8w 


ttlk 


r |(1 + w) 



16 
where 



H = !— 1 — i 1 — | \ , ep 

p 2 + ■ ■ ■ (77) 



—8mK 


r 


— 8mn 


r 


\w 2 + QAw + m 2 K 2 


r 


2 w 


m 2 K\ 


r| (1 + w) 3 



TTT / mK\r\ m K \r\ 
w = W \ —e s 



The leading term is simply the mass plus a static gravitational correction, which is m K\r\ 
to lowest order in k. The term linear in p is due to purely to gravity, since it vanishes in the 
limit k — > 0. The argument of the function w must be larger than — 1/e. This translates into 
the limit mur < W(l/e) which means that if mnr is sufficiently large, there is no small p 
expansion - i. e. there is some minimum value of p below which the Hamiltonian is no longer 
real. This situation is shown in Fig. 3. 



Fig.3 

Hamiltonian as a function of momentum p > and nr in units of m. 



The small m expansion (which is the same as the large p expansion) is also easily obtained. 
For example in the region p > 0, r > with e = 1 we find 



1 9 1 —4 + Krv , 1 16 — Aurp + n 2 r 2 p 2 , 

H = 2p+ -m 2 + ^-m 4 + - c —m 6 + ■■■ . (7i 

y p 16 p 3 128 p 5 
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The ^-independent terms are equivalent to those obtained for two free relativistic particles of 
equal mass in the small mass limit. We see that the effects of gravity modify the Hamiltonian 
to include interaction terms whose strength grows with increasing separation, as one might 
expect from the basic structure of two-dimensional gravity. 

The Hamiltonian ( [76|) describes the surface in (r, p, H) space of all allowed phase-space 
trajectories. Since if is a constant of the motion, a trajectory in the (r, p) plane is uniquely 
determined by setting H = H in (|76|). However there are two distinct sets of trajectories 
which correspond to the two real branches of H^-function which join smoothly onto each 
other. 

For Hq sufficiently small, the Hamiltonian is given by the principal branch Wq and reduces 
to the Newtonian limit for small k. Fig. 4 shows the Hamiltonian for small H . The darker 
surface denotes H(r,p) = H and the lighter surface is the Newtonian Hamiltonian. 



Fig.4 

H in the non-relativistic limit compared with the Newtonian Hamiltonian 



Once H becomes sufficiently large, there appears a qualitatively new set of trajectories 
which are not connected with the Newtonian Hamiltonian in small k. Fig. 5 demonstrates 
two branches of the Hamiltonian in the region —2 < r < — 1, —2 < p < — 1. The whole 
surface of H continues smoothly from one branch to the other. This structure is seen in 
Fig. 6 where the slice of H at constant p (p = 0) is drawn. 



Fig. 5 

Both branches of the Hamiltonian 



Fig. 6 

A slice of H at constant p 



For Hq in the intermediate range the constant energy surface intersects both branches 
shown in Fig. 7. The trajectory in (r, p) plane moves over both branches. It analytically 
continues from one branch to the other. Fig. 8 shows also this structure from another point 
of view. 
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Fig.7 



A comparison of the Hamiltonian with a surface of constant energy. The flat black surface 
corresponds to a value of H = 5. Note that it intersects both branches of the Hamiltonian. 



Fig. 8 

A slice of both branches of H at r = 1. The horizontal line corresponds to H = 15 



For a given initial condition the energy H of the system is fixed and the trajectory in 
(r,p) plane is given as the slice of H = H through the 2-dim surface H(r,p) in (r,p,H) 
phase space. Two characteristic plots are shown in Figs. 9-10 where the corresponding 
trajectories in the Newtonian theory are included for comparison. Under time-reversal, the 
trajectory for a given value of H is obtained by reflection in the p = axis. 



Fig.9 

Non-relativistic (Newtonian) and relativistic trajectories for H = 2.2. The undistorted oval 

is the non-relativistic trajectory. 



Fig. 10 

Non-relativistic (Newtonian) and relativistic trajectories for H = 3. The undistorted oval 

is the non-relativistic trajectory. 



One of the characteristics of the trajectories is that as H increases the trajectory becomes 
more 'S"-shaped. Suppose the particles start out at the same place (r = 0) with positive p. 
r will increase and p will slowly decrease. This continues until maximum separation with 
some positive value of p, where the velocity r — 0. After that r undergoes a rapid decrease, 
while p is still positive. At some value of r, p becomes zero and then it goes negative. The 
particles continue to be pulled together and r reaches 0, where p has its maximum negative 
value. The particles then overshoot the mark and start the reverse motion with interchanged 
positions. 
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ran 


r 


+ 2mK\r\w + 8w 


ran 


r |(1 + w) 



The main reason why the trajectories are 'S' shaped is the appearance of the p-linear 
term in the Hamiltonian. The canonical equations flSED and ([J7|) (or directly the Hamiltonian 
(|77D) leads to 

(p — terms) . (79) 

The first term on RHS comes from the p-linear term in H and r = does not correspond 
to p = 0. This relation between r and p resembles the relation in the theory with charged 
particles 

p = 7r% + eA • 

In this sense the first term on RHS of (|7"9|) can be said to be purely gravitational. 

6 The unequal mass Hamiltonian 

For unequal masses we set 



H = t-* + -(pi - P2) n Y 80 

2 2 K\r\ 



so that (BID becomes 



(Y 2 - D 2 )e 2Y = (S 2 - D 2 )e 2S (81) 



where S = ^(M 1 + M 2 ) and D — ^(Mi — M 2 ) with Mi = ^p\ + m\ - ep x and M 2 = 

^2 + m?, + ep 2 - For equal masses, S = Z, D = and (^l|) reduces to (ffT]) . Solving (|8TD for 
K in terms of S and D yields the Hamiltonian in the unequal mass case. 
To obtain the solution, consider the equation 

(y 2 - a 2 )e 2y = (x 2 - a 2 )e 2x a > . (82) 

This equation also has three solutions shown in Fig. 11 : the trivial solution y = x, the curve 
H-I-J-K-L denoted by W(x; a) and the curve S-T-U denoted by W(x; a). To our knowledge, 
discussion of the functions W(x; a) and VV(x; a) have never appeared in the literature. We 
shall refer to W as the generalized Lambert function since lim a ^o W(x; a) = W(—xe x ). In 
general W is also complex and multivalued, and when x is real, the function has two real 
branches shown in Fig. 11. The principal branch is analytic at x = and has a derivative 
singularity at x — W _1 (— 1(1 + \/l + 4a 2 )) beyond which it becomes complex. The other 

branch satisfies W < — 1(1 + yl + 4a 2 ) and joins smoothly onto the first branch. The full 
function is double valued for a < x < W _1 (— 1(1 + \/l + 4a 2 )). The third solution VV(a;; a) 
is a generalization of W(xe x ) in the region x < 0. 



Fig. 11 

A plot of solutions to eq. (82). The curve HIJKL is the Generalized Lambert function. 
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As in the equal mass case, the trivial solution Y = S again yields the unphysical Hamil- 
tonian H = e{p~\ — fiz). Since the physical domain of 5* is 5* > 0, the only physical solution 
is 

Y = W(S;D) (83) 

which leads to the Hamiltonian 

H=^ - + ^Pi-P2)-SyW[^(M 1 + M 2 ) ] f ^(M 1 -M 2 )] . (84) 

The expansion in k for this Hamiltonian is given by fl62|). 

Choosing also the center of inertia frame with pi = —p 2 = p and setting m = m^/mi-, we 
shall look at the trajectories. First, take a value for Hq just above the minimal (rest-mass) 
value of 1 + m and compare this to Newtonian theory in Fig. 12. 



Fig. 12 

Non-relativistic (Newtonian) and relativistic trajectories for H = 2.01 in the unequal mass 
case with m = 1. The undistorted oval is the non-relativistic trajectory. 



The trajectory is almost exactly the same as Newtonian theory, since it is the equal mass 
case. For larger m the separation between particles cannot get to be very large and the 
trajectory becomes more compact as shown in Fig. 13. 



Fig. 13 

Non-relativistic (Newtonian) and relativistic trajectories for Hq = 6.01 in the unequal mass 
case with m = 5. The narrow oval in the middle is the relativistic trajectory. 



The trajectories for smaller values of m are shown in Fig. 14, where the innermost line is 
the m = 0.9 case and the outermost is the m = 0.1 case. 



Fig. 14 

Relativistic trajectories for several values of m, where H = 2. 
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Finally in Fig. 15 the trajectories of different values of m both large and small are compared. 



Fig.15 

Relativistic trajectories for several values of m, where Hq = 4. 

7 Solution of the metric tensor 

To determine the Hamiltonian and derive the canonical equations of motion, we had only 



to solve the constraints (|i2[) and fll3|) of the system of the field equations fllQT) - (|17|) . In this 
section we shall solve the remaining equations to determine the metric and to confirm directly 
the consistency of Euler-Lagrange equations fllB] ) and (|i~7|) with the canonical equations 



derived from the Hamiltonian, though formal proof of the consistency was already given in 

i- 

Under the coordinate conditions (|23|) the field equations (|T(]|), (|TTD, ([TJ) and (^) become 



55) 



/tvriVo + AT| = (86) 
d 1 (±N o *' + Nti = (87) 

* = -2kttN + Ni9' . (88) 

In the following we shall carry out our calculations assuming z 2 < Z\ - the case z\ < z 2 
is completely analogous and will not be presented. The solution to (|87D is 

( A(f) 2 + (+)region 
N = Ae~^ = A<p 2 = < A(j>l (O)region (89) 

AcfPi (—) region 



where A is an integration constant and <f>± and O are given in (0), (0) and (|57|). Eq.(p6|) 
becomes 

N[ = -kA X '<P 2 . (90) 

The solution in each region is 



,+) region : 



N 1W = e [ A ^ L ^ - l} = e (A<Pl - l) (91) 
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(-) region : 

AT 1( _ } = - e l^e-S^-t—") - l} = -6 (A/> 2 _ - l) (92) 



(0) region : 

= ^ " &«- ,M ~i + f ^ + fC » (93) 

where we chose the integration constants in Nu + ) and Nit-\ to be (— e) and e, respectively. 
In general we should respectively replace these with arbitrary constants C + and C_ appear 
in (|9~ID and (|92|). However a lengthy calculation reveals that they are equal to (— e) and e, 
respectively. For simplicity we shall set C + = — e and C_ = e from the outset. In deriving 
iVi(o) we used fl63|). 

The continuity condition at x = z\ 

N l{+) ( Zl ) = N m ( Zl ) (94) 

leads to 

q, = -i + JiL-fr(fr-M.)_«jA 1 (95) 

\Afi ioA/i 2 L„ 'J ^ ' 

where ( |B"3[ ) and the relation L 2 + Mi = L are used. The continuity condition at x = z 2 

N 1{ ^(z 2 ) = N m (z 2 ) (96) 

similarly leads to 

Co = 1-a\^- L ^- M *K^z 2 ) (97) 
\M 2 L M 2 2 L 2 J 1 ; 

From the consistency of fl95D and (]97f) the constant A is determined as 

Li + L2 — ^fL\L 2 ^ 

and 

Mi - M 2 - + z 2 )L 1 L 2 , . 

C = r r 4 Lr r • (") 



Li + L 2 — ^fLiL 



2 



Now we are ready to check fl85l) . First we treat three regions (+), (0) and (-) separately, 
and next consider the matching conditions at x = z±, z 2 - 

For (+) and (-) regions it is straightforward to show that the LHS of ( jB5|) vanishes, by 
substituting the explicit solutions of n = x', ^' = — 40'/0, N and N^. For the (0) region the 
solutions of the metric and the dilaton field are 

71 = —^Lq, 7r' = 

7T = -|l = -|(4X - epi + ep 2 ) = ipi 
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00 

N = A<j> 2 , N^ = 2AM' 



The LHS of fl85[) becomes 



1 9 A 

LHS of (85) = /i + |f(^o) 2 -^(0o) 2 

1 . k LqL x L 2 
:Pi + 



2 8 h\ + L2 — ^-L\L 2 

which vanishes due to (|64"D. 

Consider next the ^-function part at x = Zi. Since 

7r = = — T {Pi I x — z i I' +P2 \ x — z 2 \'} — eX 
7r' = -~ {j5i5(a; - 21) +p 2 5(a; - z 2 )} 

* = \ {V\Z\b{x - z x ) + p 2 z 2 5(x - z 2 )} - - {pi \ x - z x \ +p 2 \ x - z 2 |'} 

we have 

[5 — function part of (11) at x — Zi] 

1 1 p 2 1 

= -pizx5{x - zi) - -iV (^i)-T= 1 8(x - zi) + -N 1 (z 1 )piS(x - z x ) 

Z Z Jp2 _|_ m 2 Z 

1 xt ^ J • eL ° Ll 1 

-pi<)(x - Zl) < + 



2 [ Li + L 2 - j5 fL 1 L 2 ^Jpj _)_ m 2 

which also vanishes due to the canonical equation (|66|). Similarly, the 5-function part of (11) 
at x = z 2 is zero. We thus conclude that fl8"5j) is satisfied exactly. 



As we investigated in the iterative method [§], for the consistency of (|38D we need to 
introduce a x-independent function f(t) into ^ 

^ = -41og I I . (100) 

Since in the system of the original equations (|lO"l) -(|T7l) all other equations except fll4l ) contain 
only spatial derivatives of f(t) does not contribute to either the Hamiltonian or to the 
equations of motion. Eq. (^) becomes 

— 4^ + f(t) + 2kttNq + A^-Ni = (101) 



We must check this equation in the three regions separately, with f(t) common to all regions. 
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For the (+) region, after substituting the solutions of (f> + , it, iV (+) and iVi( + ), C P- Q H ) yields 

■ , . L\M\ — L\M\ k , 

/+(t) = 2 LlMl ~ 2 i+(il - £) (102) 



This ensures that /+ is x-independent. Using the canonical equations (|64|), (p5|) , (pq) and 
(0) we get 



/+(*) 



re r 
2 L 



Li + L 2 — b j-L\L 2 



U 6 



Pi 



+ £ 2 e 



P2 



pi + mi 



(103) 



For the (-) region the calculation is quite analogous to the above and leads to f-(t) = f+(t). 
For the (0) region the calculation is rather lengthy and complicated, especially for 0o, which 
is expressed as 



Lq 2,L\ c 2,L 2 8 
e{M 1 - M 2 ) 



-^1 + L>2 — 



+ —X(p + Kr ^0 



8 W > \\M X 



T \ 1 / 2 



Mo 



z 2 e *' 



Substitution of the expressions of O > 0o> 4>o, ^ ^0 an d into (|101|) in the (O)-region leads 
to 



fo(t) 



. ( Ln La Lo \ .LiLo 



9-^0 



Li + L 2 — -rL\L 2 



LA 6 



L 
Pi 



+ L 2 \e 



P2 



P2 + m 2, 



(104) 



which is equivalent to ( 103 ). Hence /(£) is common in all regions 

f + (t) = j-(t) = Mt) ( 105 ) 

and the solution is self-consistent. 

Finally we shall directly check the Euler-Lagrange equations fllCf ) and ( |i~7D which under 
the coordinate conditions (E3) become 



z a - N (z a ) 



Pa 



+ N 1 (z a ) = 



Pa + 



Wo 



Pi + m z a 



2~ —2 dN ^ n 
Pa a ~ ~dz~ Pa = ' 



(106) 
(107) 
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Since N and Ni are continuous at x = zi, z 2 , we have 



N ( Zl ) = A<Pl(zx) 
N x { Zl ) = 1 
Then for particle 1, say, (|106|) is 



Lq Li 

~~ L x + L 2 - f LiL 2 ' M~i 
L Li 



1 



Zi= e 



eLoLx 



L x + L 2 — ^L X L 2 , /p| _|_ m 2 



(10? 



This is identical with the canonical equation (|66|). 

On the other hand, 8N /dx and dNi/dx are discontinuous at x = Zi,z 2 . The natural 
definition of dNo/dz x is 



dN _ 1 / 0jv o 



2 &r 



+ 



rr=zi+0 



dx 



x=z\—0 ■ 



(109) 
(110) 



and similarly 



0iVi 



1 f <9iVi 



2 &r 



+ 



dN x 



x=zi+0 



dx 



x=z\—0 ■ 



Kt A L\ , r T \ 

—A—(L + + L ) . 

8 Mi + 0J 



For particle 1, ( |107|) is 



Pi 



[LqL\L 2 



L\ + L 2 — ^L X L 2 



:nv. 



(112) 



;ii3) 



which is identical to (64) . For particle 2, (P-06 ) and (|107|) also reproduce the canonical 
equations flBSf ) and fl6T|). 

Thus the consistency of the solution has been completely proved. 



8 Test particle approximation 

As an interesting limiting case of (|6ll) let us try to get the Hamiltonian in the test-particle 
approximation. 

Setting particle 1 to be a test particle \i and particle 2 to be a static source m at the 
origin, namely, 



Zi = z, 

Z2 = 0, 



m x = 
m 2 = m, 



Pi =P, 
P2 = 0, 



Pi=P = P-^ 

P2 = 
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the defining equation (|6l|) for the Hamiltonian becomes 

If + /J 2 - H)(m + ep - H) = (Jp 2 + fi 2 - ep)m e^ H ~ e ^ zl . (114) 



Expanding if in a power series in ^Jp 2 + yU 2 and ep and taking only the linear terms we 
obtain 

H = m+^p 2 + fx 2 eTW - ep y- (e^ z \ - l) . (115) 

for the Hamiltonian in the test particle approximation. This Hamiltonian is expressed in 
terms of the metric tensor of the static source as 



H = m + \jp 2 + /i 2 N (z) — pNi(z) (116) 

where 

N = e^ = c J?L ( e W - i) . (117) 

The canonical equations are 

P 



yfp 2 + [i 1 

and 



N - Ni (118) 



Eq. (|118|) is inversely solved as 



p= " W + *> , ■ (120) 



We set the initial condition 

z = z = v at t = 0. (121) 
Then the initial momentum p(0) and the total energy H Q are 



fJLVQ 



1-v 2 



p(Q) = f : iio = m + v /j9(0) 2 + /i 2 (122) 



From ( |116| ) and (|122|) p is given by 



p= N 2 l _ Nl { M°) 2 + ^ ± N oM°) 2 + ^ - ( N o - N ? )/^ 2 } • ( 123 ) 

We can draw a trajectory in phase space, an example of which is shown in Fig. 16. The 
trajectory is again 'S' shaped due to relativistic gravitational effects. 



Fig. 16 
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Relativistic trajectory for a test particle compared to the Newtonian case. 



From fll20p and (|i23|) we get 



iV n 



N 1 + N Jl-(N 2 -N 2 )(l-v 2 ) 



(124) 



'(iV 2 - iV 2 ) 2 (l - vl) + [N, + N ^l - (JVg - N?)(l - v 2 ) 
Denoting the RHS as G(z), the solution is given by 

t 



dz 



o G(z) 



To lowest order in k we obtain 



and 



G(z) 



vl- 



Km , 



Km 



: t 2 + v t 



(125) 

(126) 
(127) 



which is the Newtonian motion of a test body in (1+1) dimensions. 

We shall add one comment on the form of the line element. In our canonical reduction 
we chose the coordinate conditions (^4|), under which the line element of space-time is 

ds 2 = -(iV 2 - Nl)dt 2 + 2N 1 dtdx + dx 2 . (128) 

For the case of a single static source with No and N\ given by ( |117|) , we find the coordinate 
transformations 



t = t - e\x\ + — log|2e^ |s| - 1| 
Km 



x 



X / urn 

Km \x 



(efM - l) 



(129) 
(130) 



which leads to the line element 15 



ds 2 



-a(x)dt + 



MX) 



-dx 



z,2 



with 



a (X) 



1 + 



Km , 



x 



In this coordinate frame the Hamiltonian for the test particle is 



H(z,p) = Ja(z)fi 2 + a(z) 2 p 2 + m 



(132) 



(133) 
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9 Correspondence with Newtonian gravity in (d + 1 
dimensions 

In this section we illustrate how a Newtonian limit generically arises in the (1 + 1) dimensional 
theory we consider. We compare this with the emergence of a Newtonian limit in id + 
1) dimensions. We shall compute the Newtonian limit (s) by considering the one graviton 
exchange potential (keeping in mind that there are no propagating gravitons in two spacetime 
dimensions) . 

We begin by extending the theory in (|I|) to d + 1 = n dimensions and coupling N scalar 
fields, which yields 

L = {*R + ifTV^V,*} -\j2V^9 {g" u <f> a ,^a,u + mlrt) (134) 

where k 2 = 32ttG. Defining the graviton field and the dilaton field ip via 

SVu/ = V)iw + Ktl^ u * = 1 + /C0 (135) 

gives 

L = ~{d x h^d x h^-d x h^d x h\-2d,h^d x h Xu + 2d,h^d u h x x } 

-2(crV- d,hVd^ + d^d^-~Y,(W0a,, + rnlti) ■ (136) 



for the free Lagrangian density following from ( |134| ). Redefining the dilaton field 



QiXQV 
□ 



i> = V + h% - — V (137) 



allows us to express (|136|) as 

L = ~{d x h^d x h^-d x h^d x h\-2d^d x h Xv + 2d,h^d u h x x } 

- (d" V - W (d x h\ - d»h\) + d^d^ -W Uf<p a , „ + mlfi) .(138) 

The field ip decouples from the Lagrangian and we shall not consider it further. 

The free Lagrangian density of the graviton is obtained by simplifying the first two terms 
above 

L 0g = ~\ {d x h^d% v + d x h^d x h\ - 2d^d v h\] + d v h^B» + X -B^ (139) 

where we added gauge fixing terms in the form of a Lagrange multiplier field B^. 
Eliminating B^ from its field equation leaves us with the Lagrangian 

Log = ~\ {d x h^d x h flu + d x h^d x h% - 2d^d„h\} - d»h^d x h» x . (140) 
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whose canonical quantization we shall now undertake. 
Temporarily setting the scalar fields to zero, we obtain 

□V + V^h\ - ri^d x d p h x P - d^d v h\ + d p d x h vX + d v d% x = . 

for the graviton field equation. Its trace is 

n — 2 

Uh\ = d x d p h xp 



n 



implying that ( |141j ) becomes 

□V - - W\ + ^ + = o . 

Taking the d u derivative of (|143|) leads to 

□0" v = o . 

This, along with the D'Alembertians of (|143|) and (|144[) respectively lead to 

□ 2 V - d p d u ah x x = and D 2 /i A A = 

which finally implies 



□ d V = 



a relation characteristic of n > 3 Lagrangians. 
The conjugate momentum is 



7T 



doh^+rT (d h x x + d A /i A )+r^° - d x h vXS j+rf° Q<9^ A A - d x h p) 

which implies 

<9o^-oo 



2 2 

~2 n ° l + 4 dih °° ~ 4 9lhlj + 2 9jkij 



dnh 



O'Hj 



n 



n 



The equal-time commutation relations are 



{h p „(x), h Xp (y)} 



cq 



vr^(x), 7T Xp (y) 



. 



<■</ 



implying that the commutators between hu V and 9o/iAp become 



i cq 



(141) 
(142) 

(143) 

(144) 

(145) 
(146) 

(147) 



(148) 



(149) 



(150) 
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where 

%v = Vnv + VfiaVvo ■ (151) 

The proof of this relation is given in the Appendix. 
The solution to (|146|) is 

K v {x) = -J d n -h - z)B z h IIU (z) - J d n - l z D^(x - z)d z Q Dh pu (z) 

- I d n - x z D (x- z)d z a 2 h^ u (z) . (152) 



where D^ n \ Z)W and D are defined via 



D {n \x) = "T^pi / d n k e(k )S(k 2 y k * 



D {n \x) = —j^t—[ I d n k e(k )8'(k 2 )e ik3 



(27T 

~ (n) j /■ 

d (x) = -^pr y <™ e(*b)<5"(*V 



We next need to express all of h^, doh^, ndoh pu , d 2 /?.^ and D 2 doh fJiU in terms of the canoni- 
cal variables and calculate commutators at equal-time. This rather lengthy and complicated 
calculation is given in the Appendix. 

From (|152|) and the equal-time commutators, the commutator among the components of 
Ufa, at two arbitrary space-time points can be calculated : 

[h^(x), h Xp (y)} 

+ ^ {-Tlpxdudp ~ V^dudx - r]uxdpd p - r] up d^d x + - (r]pud x d p + i] Xp dpd u )^ D {n) (x - y) 

if 2 \ ~ («) 

+ - 1-- )d p d u d x d p D (x-y) . (153) 



n 

This expression is valid even when n = 2. The graviton propagator is 

(0|T(V(*)M2/))|0> = J d n k e *(»-»)^^ (154) 

where 

2 



1 r 4 i 

2^2 | — VfiX^vkp — f]ppk u k x — Tj uX k p kp — f] up k p k x + — \Tjp, v k\k p + f] Xp kpk v ) > 



+ ,l-§)4£*£. (1B6) 
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We turn now to the scalar fields, whose Lagrangian density to lowest order in the graviton 
coupling is 

L int = ~\ {\v" u O' a 0,a + m 2 <f) - <w} V • (156) 
The one graviton exchange diagram yields the S'-matrix element 



_ AmG n , Q Q Q x-i/2 
~ (2tt)™- 2 \ PlP2qiq2 J 



P2Q2 ~ ^V Xp (P2 ■ Q2 + mj) 



x ^P 5 (n) (pi+p2 _ qi _ q2) 



where p^, q£ and /c M are the four- momenta of the initial particles, the final particles and 
the transferred graviton, respectively. This result is valid for n = 2 also. In the lowest order 
and the static approximation T-matrix element is 

Tn = -4 (l - -) 7^4^? • (157) 

whose associated potential is V = J d n ke~ tkx T(k) in n dimensions. 
The T-matrix elements for n = 2, 3 and 4 are 

2G 2 mim 2 

T 2 = — (158) 

8 G 3 mi m 2 
T3 ~ "3 (>) ~k^~ (159) 
3G 4 mim 2 

T4 " " (2^ ' • (160) 

and the corresponding potentials are 



V 2 = 27rG' 2 m 1 m 2 r (161) 
V 3 = 2 ^G 3 ) mim 2 log r (162) 

K 4 = (163) 

By identifying the gravitational constants as 

4 3 
Gn,2 = G 2 Gjv,3 = -C 3 Giv,4 = -G4 (164) 

we get the correct Newtonian potentials in each dimension. 

The above results are in strong contrast with (d+1) dimensional general relativity, whose 
free Lagrangian density is 

L 0g = -\{d x h^d x h^-d x h%d x h%-2d^d x h Xv + 2d^ v d„h\} 
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where gauge fixing terms have been added. A computation analogous to the one above gives 



(o|T(V(z)My)|o) 



where 



X^xp = rjpxrivp + rj^rjux - 



V 

jn^ c ik{x-y) ^liu,X P 

k 2 — ie 



n-2 



The S'-matrix element of one graviton exchange diagram is 



S 



AixiG n 
(2tt)™- 2 



-1/2 



Q /3 

P2Q2 



~T] aP (p2 ■ q2 + m 2 2 ) 



x„ 



which in turn yields the T-matrix element 

T , 4a 



Q2 ) 



n 



3 m\m,2 



(2tt)"- 



/? 



A' 2 



(165) 



in the static appoximation in n dimensions. 
The potential for n = 4 is 

= 



(166) 



in agreement with ( |163j ). However the potential for n = 3 vanishes, and the potential for 
n = 2 diverges. This latter situation can be dealt with by setting G n = (1 — | )G% and taking 
the n 



2 limit |13| , which yields the 2-dimensional T-matix element 

2G2^1^2 



T 



k' 2 



and potential 



= 2nG2Tn\m2r 



(167) 



(168) 



Unlike general relativity, 3-dimensional dilaton gravity includes the Newtonian potential 
in any dimension, once the gravitational constant is appropriately rescaled. In this sense the 
theory of gravity (|TJ) we consider is a relativistic extension of Newtonian gravity in (d + 1) 
dimensions. General relativity, on the other hand, does not include Newtonian gravity in 
(2 + 1) dimensions and is empty in (1 + 1) dimensions. In the latter case an appropriate 
rescaling of Newton's constant yields the theory (IT]) in the n — > 2 limit [13]. 



10 Conclusions 

We have obtained an exact self-consistent solution to the 2-body problem in a (1 + 1) 
dimensional theory of gravity with a Newtonian limit. To our knowledge, this is the only 
exact relativistic 2-body solution of this type. We are able to explore all possible limits of 
this solution, including large and small gravitational coupling and/or mass and/or momenta. 
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A natural extension of what we have done would be to attempt to solve the N body 
problem. It would also be of interest to couple other matter fields (e.g. electromagnetism), 
and to investigate the extent to which our methods are applicable to other dilaton theories 
of gravity. 

Finally, and what is perhaps most interesting, is to quantize the degrees of freedom of 
the 2 body system we consider based on the Hamiltonian given in fl6ip. The quantum theory 
based on (refeqHl) is a quantum theory of gravity coupled to matter whose slow-motion weak 
field limits should be straightforwardly comparble to that of the non-relativistic mechanics 
of 2 particles in a linear confining potential. As such it should offer interesting insights into 
the behaviour of quantum gravity. 
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Appendix: Commutation relations 

We use the notation 




A 



A 



dxh^ 



) 



J eg 




(169) 



from which we shall now deduce various commutators of interest. 



Taking the trace of ( |169| ) implies the relations 



[V: nd ° hQ * ~( n ~ = ir)^5 {n - l \x - y) 

2 % 
[hp,, d h a a - d a h a ] eg = - [hftv, d h 00 } eq + -ri^ u 5 {n ~ l) {x - y) 

[V> d o ha a\ eq = ~ — - tV> do^ooL + - , n^ {n ' 1 \x - y) 



(170) 



from which it follows that 




(171) 
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The (k, a) = (0, a) component gives 

[h,a„ d h a} eq = l - (Vvoriva + VfMvVvo) 5 (n_1) (a; - y) (172) 
which yields in turn the commutation relation 

K<j\ f, n 

= \{~2 ^ Kr]va + VvvVvk) + \ {VpkVv* + V^Vuk) - -VuvfjK*} 5 {n ~ 1] (x - y) (173) 

We now consider the problem of expressing Oh^, nd h^ u , D 2 h^ u and D 2 d h^ u in terms of 
h/j, u and d^h^ (or equivalently the canonical variables). The components of the first-order 
field equations ( |141| ) give 

□/loo = ^Ahu - ^didjhij 



ah 0i = ^didohjj - ^djdohij - ^didjh 0j + ^Ah 0i 
4 

Dhij = --Sijdkdoh ok + did h j0 + djd h i0 

+-5ij (3d k dih k i + 2Ah 00 - Ah kk ) 
n 

—didjhoo + didjhkk — didkhjk — djd k hi k 

and 

^o^oo = ^didjhij + A/ioo - -Aha 

doh Qi = —dido hjj + - dj d + dj h 0j + ^Ah 0i 

2 ^ 
9 Q hij = —5ijdkdoh Qk — did hj — djd h i0 

--S i:j (3d k dih k i + 2Ah 00 - Ah kk ) 
n 

+didjh o - didjh kk + did k h jk + djd k h ik + Ahij 
From flT7|) and (|T75l) we get 

n d h 00 = ^Adohu - ^didjdohij 



(174) 



(175) 



ad h 0i = Ad h 0i + (l - didjd h 



U oj + (~ ~ l) Ad i h oo + (l - - ) A0,h n 



-Adjhij + f~ - 1 ) 0,0,01,11,1, 



(176) 
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Ddohij = -Sij (Ad h kk + d k did h k i + 2A<9 /i o - AAd k h ok ) 



n 



1 1 

-didjdoho - ^did k d h jk - -djd k d h ik + didjd k h ok 



whereas ( 145|) implies 



O 2 h 00 



u 2 h 



0/ 



and 



a 2 hij 

a 2 d h 00 
n 2 d h 0i 
a 2 d hij 



1 - -) {-2A + A 2 h 00 - l -A 2 h d + -AdidjhiA 
1 - -) di |a<9 ^oo + g Adohjj - 2Ad j h 0j + -djd k d h jk \ 
1 - - J didj <,-2d k d h 0k + Ah 00 - -Ah kk + -d k dih H > 



1 - ^) | A 2 d h 00 + ^A 2 d h u + -didjdohij - 2A 2 d l h. 
1 ~ D 9l {~ 2Ad A h oj + A 2 /i 00 - \A 2 h ]3 + ^Ac^*} 
1 - ^ |a9 /i o + ^A^/ifcfc + ^d k did h kl - 2Ad k h 0k 



(177) 



(178; 



To calculate the commutator [h^x), h\ p (y)] at two arbitrary space-time points, we first 
express h Xp (y) as 



h\ P {y) 



- x z D^\y - z)d z h Xp (z) - J d n - l z m{y - z)d z nh Xp (z) 



-\d n - x zD (y-z)B z n 2 h Xp (z) 



(179) 



\h pu (x), h Xp {y)\ 

= -J d n ~ x z : D^(y - z) [h pu (x), d h Xp (z)] eq - d D^ n \y - z) [h^x), h Xp (z)] eq 

+ D^(y-z)[h pu (x), d nh Xp (z)] eq -d D^(y-z)[h pu (x), Uh Xp ) eq 



~ (n) r -, — (n) r 

+ D (y-z) h pu (x), d a 2 h Xp (z) - d D (y - z) h pu , a 2 h Xp (z 



eq 

(180) 



We shall refer to the 6 terms on the RHS as term 1, term 2 • • • term 6 respectively. From 
the commutator (|150|) , terms 1 and 2 are 



(term 1) 



jv^p + V^pVuX + VpWup + rj pp rj uX - -77^77^ j D (n \x - y) 



(term 2) = . 



(181) 
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For terms 3 and 4 we shall calculate the equal-time commutators [h pu (x), d nh\ p ] and 
[hfj, u (x), nh\ p (z)] eq by using the expressions ( |1 74[ ) and ( |176| ). Putting terms 3 and 4 together 
gives 

(term 3 + term 4) 

i 

= j {-Vnxdvdp - fj^pdydx - fjuxd^dp - f] vp d p dx 
4 

+ - {VfMudxdp + r]x p dp,d v ) - 2r] ll0 r]uodxd p - 2r]xor] pQ d p d v 

1 - (vfjoTjxodvdp + rinoVpodudx + VuoVxod^Bp + r] u0 r] p0 d p dx) j D (n) (x - y) 
i f 4 

+ 4 U ^oVvoVxp + VxoV P o%u) 

- (v^oVxaVvp + VnoVpoVvX + VuoVxoVpp + VmVpoVpx)} D {n) (x - y) (182) 
where we used the relation 

AD {n) = dlD (n) + D {n) . 

Similarly by calculating the equal-time commutators [h pu , do0 2 hx p } eq and [h pu (x), n 2 hx p (z)] eq 
(making use of the expressions ( |177| ) and ( |178|) ), we have 

(term 5 + term 6) 

if 2 \ ~ ( n ) 

= - ( 1 - - J d p d u d x d p D (x - y) 

+7T (l - - ) {VnoVxodvdp + i] fl0 r]pod,ydx + r] u0 i]xodpd p + i] M r] p0 d p dx 



2 V n 

where we used 



2 V n 

+ VpoV»odxdp + VxoVpodpdv 

Vp.oVv 

+ ^Vp,oVuoVxoVpodo} D {n) (x - y) 
Ul--) ^oVu ^poD {n \x - y) (183) 



~ (n) ~ (n) _ , . 

AD = d 2 D + 



Summing up all the terms we get 
[h pu (x), hx P (y)] 

= l ~ [p^xVup + VppVuX ~ -VpuVxp^j D (n) (x - y) 

+^ {-Vnxdudp - r] pp d u d x - r] uX d p d p - r] up d p d x + ^ (v^dxdp + T] Xp d p d u )^ D (n) (x - y) 

if 2 \ ~ ( n ) 

+ 2 (l - -J d p d u dxd p D (x-y) . (184) 
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